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Abstract 

We discuss random matrix models for the spontaneous breaking of both chi- 

ral and color symmetries at zero chemical potential and finite temperature. 
> 
J^ ' Exploring different Lorentz and gauge symmetric color structures of the ran- 

f^ ' dom matrix interactions, we find that spontaneous chiral symmetry breaking 

^H ! 

On ' is always thermodynamically preferred over diquark condensation. Stable di- 

ON 
r~| ' quark condensates appear only as SU{2) rotated chiral condensates, which do 

Q^' not represent an independent thermodynamic phase. Our analysis is based on 

1-^ ' general symmetry arguments and hence suggests that no stable and indepen- 



ON 
ON 



o 



X 



dent diquark phase can form in QCD with two fiavors at zero quark chemical 



j^ ■ potential. 



I. INTRODUCTION 

The exploration of the phase diagram of QCD is of fundamental interest for understand- 
ing the nature of strongly interacting matter at high temperatures and densities and its 
implications for the physics of heavy-ion collisions and neutron stars. One central aim of 
studies of the phase diagram is to determine how its global structure is shaped by underlying 
symmetries, their spontaneous breaking, and their interplay. 

The current picture of the QCD phase transitions is shaped by various symmetry breaking 
patterns [0,0. For QCD with two light flavors, lattice simulations |^ indicate the restoration 
of chiral symmetry for temperatures above T ~ 160 MeV as one moves along the zero baryon 
density axis. This phase transition is of second order. Along the zero temperature axis, one 
expects chiral restoration to proceed at high densities through a first-order transition [Q. Be- 
cause QCD is asymptotically free, it is also likely that a transition to a plasma of deconfined 
quarks and gluons occurs at the same high densities required for the screening of color inter- 
actions. The ground state of this weakly coupled plasma may be dominated by quark-quark 
correlations as suggested by early and recent works @,|,^|[: the attractive quark-quark in- 
teraction may lead to the formation of quark Cooper pairs, which spontaneously break color 
gauge symmetry. Nambu-Jona-Lasinio models [^,0 and instanton-based calculations |]^,|^ 
suggest sizeable diquark condensates with {qq) ~ 100 MeV at quark chemical potentials 
/i ~ 300 MeV or higher. 

Current microscopic models of color superconductivity |@-§ are all implemented at mean- 
field level and, although different in their treatment of the primary interaction among quarks, 
have the following schematic phase diagram. Diquark condensates appear at the low densities 
usually associated with broken chiral symmetry, but they are unstable. They are found to 
be stable in the high density, chiral symmetric phase. Between these two limits, a state 
where both chiral and color symmetries are broken appears either as a saddle point of the 
thermodynamical potential , or as a global minimum [^ . Given the range of treatments 
of the QCD interactions, it is desirable to understand the extent to which predictions of 



diquark condensation are model- independent. This leads to the question of determining 
how the phase diagram is shaped by global symmetries on one hand and by the detailed 
dynamics of the interactions on the other. 

We wish to address this question by means of random matrix methods applied to SU{3) 
QCD with two flavors. As in lattice calculations, significant complications arise from the non- 
Hermitean character of the Dirac operator in models with a finite quark chemical potential. 
Because our main goal is to set the stage for studying diquark condensation, we shall confine 
our present exploration to Hermitean matrix models with zero chemical potential. 

Before introducing our method in detail, let us briefly recall the principles and moti- 
vations of existing random matrix models for chiral symmetry breaking since our for- 
mulation will follow in a similar spirit. Chiral random matrix models (xRMM) study the 
spectrum of the Euclidean Dirac operator in a theory that respects all global symmetries 
of QCD, but in which the detailed dynamics of the interaction among quarks is replaced 
by averages over random matrices. The chiral order parameter, < tptp >, can be obtained 
from the density of the smallest eigenvalues of the Dirac operator with the aid of the Banks- 
Casher formula []TU|, and properties of the phase diagram can thus be determined. The 



primary advantage of such an approach lies in its universality |]rT|. The statistical prop- 
erties of the small eigenvalues, and thus mean-field predictions regarding the chiral phase 
transition, depend only on the symmetries of the Dirac operator. They do not depend on 
the particular procedure used to average over the random interactions. In particular, chiral 
random matrix models reproduce the mean-field exponents expected for QCD with two light 



flavors 12 



In this paper, we consider a random matrix model in which both chiral and diquark 
condensation can take place and study the competition between these two forms of order. 
The paper is organized as follows. In Section ^, we start by constructing a class of in- 
teractions that includes chiral and color symmetries. The dynamics of the interactions are 
described by Gaussian averages over random matrices. In Section |IT|, we study the four- 
fermion potentials produced by integration over all possible choices of Lorentz invariant and 



gauge symmetric interactions. Particular combinations of these choices give rise to four- 
fermion interactions of the form encountered in instanton and NJL models. In Section 0, 
we introduce auxiliary variables in the usual way, derive the thermodynamical potential, 
and determine the resulting phase diagram in parameter space. As we discuss in Section |V], 
strong dynamical constraints on the coupling parameters prevent the exploration of certain 
regions of the phase diagram. It results that the spontaneous breaking of chiral symmetry 
is thermodynamically preferable to diquark condensation in all cases. Stable diquark con- 
densates can appear at most as SU{2) rotated chiral condensates, which belong to a broken 
chiral symmetry phase. This case is considered in detail in section |VII| . Our analysis of 
the phase diagram clearly distinguishes between features that directly follow from global 



symmetries and those that are model-dependent. We illustrate this point in Section y7\\ , 
where we compare our phase diagram to the predictions of NJL-model studies. 

Not all the interactions to be considered here are directly related to QCD. However, 
our arguments are based on simple symmetry considerations that also apply to QCD. This 
suggests that no mean-field treatment of QCD with two light flavors can find stable diquark 
condensates at zero chemical potential. 

II. CONSTRUCTION OF A RANDOM MATRIX MODEL 



We now turn to the construction of a random matrix model suitable for the exploration 
of diquark condensation as well as chiral symmetry breaking. Since the order parameter 
for color superconductivity is overall antisymmetric in flavor, spin, and color, these quan- 
tum numbers should explicitly be included in the model. The extended block structure of 
the resulting random matrix model is somewhat complicated. This is an unavoidable conse- 
quence of our desire to permit competition between two distinct mechanisms for spontaneous 
symmetry breaking. Fortunately, the description of the QCD partition function that will 
emerge from our analysis is very simple. While we will discuss in later sections how our 
results naturally follow from global symmetries, we now focus on the detailed description of 



the model. 

A natural starting point is to extend the xRMM by decomposing each chiral subblock 
of the single quark propagator into matrices with a flavor, spin and color block structure. 
This results in a vector theory which, as we show later, does not support diquark conden- 
sation at zero quark chemical potential. We wish to explore further possibilities of matrix 
models that build correlations in the diquark channel, and study the resulting symmetry 
breaking patterns. We are forced in the process to give up some of the symmetries of QCD. 
However, as we discuss later in length, our exploration brings insight into the question of 
how global symmetries shape the phase diagram. 

In the following, we only consider models for which the single quark propagator is Her- 
mitean. This choice ensures that the chiral and diquark order parameters can unambiguously 
be related to the spectral properties of the Dirac operator. We further impose the Dirac op- 
erator to have a 4 x 4 Lorentz invariant block structure in the vacuum, and a color symmetric 
Nc X Nc block structure. Once the flavor, spin, and color structure have been specifled, one 
is left with N x N block matrices. It should be stressed that the physical meaning of the 
matrix size A^ here is different from that in xRMM. Because we explicitly include color 
degrees of freedom, we can no longer make analogies with instanton overlap matrices, nor 
can we directly relate A^ to the density of zero modes. Ultimately, we will discuss various 
interactions by comparing the four-point potentials induced by the random background in 
the chiral and diquark channels. Hence, we regard the thermodynamic limit A^ ^ oo as a 
way to obtain mean fleld results from these potentials. The connection to zero modes is not 
lost, however, as one still expects that order parameters depend only on the dynamics of 
the smallest eigenvalues of the Dirac operator. 

Assuming a zero vacuum angle, the flnite temperature partition function has the form 
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where 1 and 2 denote flavors and where ipi, ^|, tpj, and ■ip2 are independent Grassmann 
variables which are not related by complex conjugation or transposition. 7i is intended to 
mimic the interaction of a single quark with a gluon background, and VH is the random 
matrix measure, which we will write explicitly below. We have adopted a Nambu-Gorkov 
representation for the flavor structure by transposing the single quark operator for flavor 2. 
This representation enables the random matrix interactions to build correlations between 
states with the same baryon number. We also include external mass and diquark parameters 
m and t]. These parameters permit the selection of a particular direction for chiral and 
color symmetry breaking and are to be taken to zero in an appropriate order at the end of 
the calculations. The mass m is a color diagonal matrix, associated with the chiral order 
parameters (V'lV'i) ^^id {11^2 '^2)- The complex quantity t] is to be associated with the diquark 
order parameters (V'J-PaV'i) ^ind ('?/'iPa^2)- ^^ consider here diquark condensates in a 3 
state, i.e., in a scalar state which is antisymmetric in spin and in the condensing colors. 
Hence, Pa = C'')^{i\2)-, where C is the charge conjugation matrix, and 1X2 is antisymmetric 
in colors 1 and 2. 



The temperature T enters the model in the usual manner [Q through the chiral block 
matrix 



r 



' ttT^ 



(2) 



where only the lowest Matsubara frequency has been retained. In Eq. (|1]), T appears twice 
with a positive sign. This corresponds to taking opposite Matsubara frequencies for each 
flavor. This choice is justified for s-wave pairing which, in microscopic models, leads to a 
homogeneous order parameter, < 7/'J(a;)PA^i(0) >~ A, and thus mixes Fourier components 
with opposite four-momenta. 

We write the Dirac operator Ti as an expansion into a direct product of the 16 Dirac 
matrices T^, times A*"^ color matrices: 

16 N^ 

nMak;.jPl = E i^c)^ E Kp (^a:),^ • (3) 

C=l ■' a=l 
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Here, the indices {X,k), {i,j), and (a,/3) respectively denote chiral, spin, and color quantum 
numbers. The indices {k, I) run from 1 to A^. The terms A'* represent the color matrices A" 
when a < N^ — 1 and the diagonal matrix {6c)ai3 = ^afs when a = N^. The normalization for 
color matrices is Tr[A"A^] = 26ab and Tr[5^] = Nc] the normalization of the Dirac matrices is 
Tr[rcrc/] = 4:6cc'- In the remainder of the paper we will explicitly keep track of N^ factors, 
but we will eventually be interested in the case Nc = 3. 

The random matrices A^^ represent the fields mediating the interaction between quarks 
and are thus real. Their detailed character depends on the chiral structure of Tc- In a chiral 
basis, the Dirac operator for vector or axial interactions has the form 

I W^ 



n 



V,A 



yiyt y 



(4) 



where, according to Eq. (g),W= (rv,A)RL Ea A"^r£ and W^ = (rv,A)LR Ea AM^"^. Since 
Ti is Hermitean, we have Alr = (^rl)"^- I^^ the cases of scalar, pseudoscalar, and tensor 
interactions, the Dirac operator takes the form 



n 



S,P,T 



yO Yj 



(5) 



where X = {Ts,p,t)rr EaA^A"^" and Y = {Ts,p,t)ll Ea^"^^" are Hermitean. For the 
scalar and pseudoscalar cases and for some of the six tensor matrices, X and Y are spin- 
diagonal and thus must be described by real symmetric matrices A'-''"'. As a consequence of 
rotational invariance, the matrices A"^"- associated with spin off-diagonal components of the 
tensor interaction must also be real symmetric. 

The vector interactions mimic single-gluon exchange and respect the global symmetries 
of QCD. All other interactions break at least one of them. Axial interactions are absent in 
QCD and, in a four-dimensional field theory, may spontaneously break parity. Scalar, pseu- 
doscalar, and tensor interactions, break axial symmetry explicitly and may break isovector 
symmetry spontaneously through the formation of {4'ji''Pi,) and {iPl'4'r) condensates. This 
latter condensation can actually be avoided by retaining a finite mass parameter m until the 



end of the calculation so that it is always therniodynaniically preferable to break symmetry 
in the axial rather than vector channel. We could extend our exploration to interactions 
containing the generators Tj of SU{2) flavor and study the ensuing effective four-fermion 
potentials. Here, however, we will not adopt such an approach. It will be seen from the dis- 
cussion below that the addition of further flavor symmetry will not change our conclusions. 
Having specified the random matrices A^^, we now define the measure VH 



I Ca Xk, ) L Ca \k 



(6) 



where T>A'^^ are Haar measures. The index j3c is jSc = I for real matrices (C = V, A) and 
/3c = 1/2 for real symmetric matrices [C = S, P, T). In order to respect color symmetry, the 



members of each A*"^ — 1 color multiplet share the same variance. In other words, S 



Ca 



^c\o 



for 1 < a < N^ — 1, and T^ca = ^c\s for the singlet component (a = A*"^). Further, once the 
color channel a = S,0 has been chosen, we respect Lorentz invariance in the vacuum by 
requiring a common variance T,t\s,o for each of the six components of the tensor interaction. 
The four components of the vector and axial interactions are characterized by the variances 

^v\s,o and Sa|5,o- 

III. TOWARDS A NON-LINEAR SIGMA MODEL 



We can solve the model of Eq. (^ exactly by standard methods. The first step is to 
perform the Gaussian integration over the matrix elements. This leads to a four-fermion 
interaction y and a partition function of the form 



Z{T) = I V^lV^iV^l Vi)l exp 



y + i 



' 7/;j > ' 



V. 
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(7) 



By introducing auxiliary variables, one can then derive a non-linear sigma model, the saddle 
point of which provides the exact solution to the original model in the limit A^ ^ oo. In 
this section, we concentrate on the fate of the diquark and the chiral condensates and thus 
consider only the relevant auxiliary variables. While it is largely a technical matter, it is 
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important to note that the four-point interaction has a shghtly different structure depending 
on whether the interaction is either vector or axial or is either scalar, pseudoscalar, or tensor. 
We now discuss these two cases separately. 

A. Vector and axial interactions (C = V, A) 
Integrating out the real matrices A^f^ produces a four-point interaction of the formQ 

3^c = -E^^E(4T)^ (8) 

where the quark current J^"' is 



^^ 8NT,(j^ ^,^ 



- i^im in^n^r r^u - 4iu (rL)^(A^)^ v^^v (9) 

and where Y^ = 7^^ and F^ = ^7^7^ for vector and axial interactions, respectively. When 
taking the square of Jj^l^ in Eq. (^, we need retain only those cross-terms relevant to chiral 
and diquark condensates. This becomes clear after suitable Fierz transforms of the Lorentz 
and color operators. Consider, for example, the four-point interaction resulting from the 
A^^ — 1 color multiplet. A'* = A". Denoting the Fierz coefficient in the chiral channel by f^, 
the cross-term between the first and second terms in Eq. (^ gives 

E E (ri^x)., Kp X (F^^),^., \l,^, = f^ 5,,.5^p, X 5,,5^.p + ■ ■ • , (10) 

/i=l a=l 

and leads to a four-point interaction of the form 

y- + jfkr- E ^Rki^iRk X E ^Li^iLi , (11) 

C\0 k I 

which is relevant for the formation of a chiral condensate < ip'^ip >. Similarly, the cross-term 
between the first and fourth terms in Eq. (^ can be projected onto a product of diquark 
bilinears by means of the transformation 



^ Note that the product Tr[^^2^LR] appears twice in the measure; hence the factor 8. 



4 iV2-l 
^l=l a=l 

where /^ is now the Fierz coefficient in the diquark channel, while e = (C7^)j^j^ = (C7^)ll 
and [1X2)013 are respectively spin and color antisymmetric tensors. This transformation leads 
to the four-fermion interaction 

C|0 fc I 

which is relevant for the formation of a diquark condensate < t/'JPaV'i >• 

Transposing the Dirac operators in Eqs. ( pUj ) and (|T^) provides us with the relations 
needed to transform the cross-terms of the third and fourth terms and of the second and 
third terms in Eq. (P). Taking into account the contributions from the N^ — 1 color multiplet 
and the color singlet leads then to the total four-point interaction 
1 1 

3^C = +TT7 E ^^2— E {/x (^Ififc^liJfc X ^lu^lLl + '4^2Rk^*2Rk X ^L^2L«) 
^^^ a=0,S ^C\a k,l 

+ /a (V^Lfc (^^^2) i^*2Rk X ^L (^^^2) ^iLl + ^lu (^^^2) ^*2Ll X ^2Rk (^^^2) V^i/Jfc ) } • (14) 



B. Scalar, pseudoscalar, and tensor interactions (C = S, P, T) 

The random matrices representing the interaction in the scalar, pseudoscalar and tensor 
channels are real symmetric. Their integration produces the four-fermion interaction 

3^c = - E T^^ E (4T + Jik'? , (15) 

where a again denotes the color channel [a = O, S) and, in the case of a tensor interaction, 
jj, = 1, . . . , 6 labels the six antisymmetric Dirac matrices. The quark current here is 

Jkt = 4Rk (r&)RRA'^ i^lRl - i^2Rl (r&)RR(A'^)^ r2Rk + {R - L} . (16) 



Once again, only some of the cross-terms in Eq. (|T^) contribute to the chiral and diquark 
channels. After the necessary Fierz transforms, these terms lead to the four-fermion inter- 
action 
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a ^^^^C\a k,l 

+2 fl 4^, {teX^) r2Rk X i^Ini (^^A^) i^^^, + {R ^ L}} . (17) 



IV. PHASE DIAGRAM IN PARAMETER SPACE 

We are now in a position to derive a non-linear sigma model from the total four-point 
interaction, y = I]c3^c with yc given by Eqs. (|14D and (0). We first write 3^ as the sum 
of squares of fermion bilinears and then transform the quartic terms into fermion bilinears 
by means of the Hubbard- Stratonovich formula 

r 0-2 

exp(-AQ2) ^ / daexp{-—+iQa), (18) 

which introduces an auxiliary field a. A number of simplifications now follow if one is 
interested only in those correlations among quark states which contribute to chiral and 
diquark condensates. We anticipate that the saddle points of our non-linear sigma model 
will be described in the space of auxiliary variables by dynamical masses that are real and 
independent of chirality and flavor and by complex pairing potentials that are independent of 
chirality. In the following, we restrict our attention to such choices and introduce accordingly 
real chiral and complex diquark fields a and A. 

Before proceeding, we also consider the possibility that a spontaneous breaking of the 
color symmetry affects the chiral order parameters. This possibility clearly arises in the in- 
stanton model of color superconductivity of Carter and Diakonov [|^. These authors pointed 
out that the masses of the quarks with the condensing colors must be different from those 
with the transverse colors if the Dyson-Gorkov equations are to close. To take this effect 
into account in the present formulation, we include the projections of the original interaction 
onto a chiral channel with color quantum numbers described by the diagonal generator As. 
This introduces an additional Fierz constant /g, a term /§ 5jj'(A8)a/3' x Si>j{Xs)a'f3 on the 
right side of Eq. (0), and a chiral field ag- With these additions, the chiral fields couple to 
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quarks in the color diagonal combination ip'^^cr + erg As)'?/'. This can also be written as ip'^ac'ip 
where cXc = diag((Ti, Ui, a^) with 



CTg (Tg 

(Ti = cr H — -= and (73 = a — 2^^. 
v3 v3 



(19) 



This final form of the fermion bilinear makes the separation in the chiral fields explicit. We 
follow this pattern in the external masses by taking m = diag(mi, irti^m^). 

The derivation of the sigma model potential is now straightforward. Keeping only the 
three auxiliary variables o"i, o"3, and A, the total four-fermion interaction y becomes 



+ 



^^J.^ < 



\i^2R/ 



^ ^i^J 



V 



/ 



V ^*2R J 



+ {R^L} 



(20) 



exp 3^ ~ / da^dasdA exp [-4A^ (^| Ap + B{2ai + (73)79 + C{ai - a^^f /^) 

-Oc — A£:A2 
A*£:A2 (Jc 

where we have explicitly factored out chiral and spin factors in the terms quadratic in 
auxiliary variables. The coupling constants A, 5, and C are given as 



^-2 E 



S2 



-1 



5-2 E 



Ca ^Ca 



a \ -1 / fa \ -1 

^-2 E^^ 



.Ca ^Ca, 



(21) 



.Ca ^Caf 

Inserting then Eq. (pO| ) in Eq. (|^, we obtain the partition function of Eq. (|l]). Integrating 
out the fermion fields, we have 



Z{T) ~ f daidcxsdA exp [-4iVfi((7i, (73, A)] 



(22) 



where the thermodynamic potential Q is 



fi(ai, a3, A) = A\A\^ + B [iS^al + (32CTia3 + (3sal) - {N, - 2) log[(a3 + msf + rr^T^] 

-2 log[(ai + mi)2 + I A + r/p + vr^T^], (23) 



and where 



A^^ + ic/fi, p,^^-^C/B, /33^1 + lc/i?. 



(24) 



Again, the external masses mi and m^ and the parameter rj are to be taken to zero at the 
end of the calculation. 



12 



The thermodynamic stable phases are given by the minima of the potential fl in Eq. (p3D. 
In practice, one does not expect any condensation in repulsive channels; the global minimum 
must give zero for the corresponding variables. In the following we therefore consider the Q 
as given in Eq. ( p3D only in cases where the chiral-As channel is attractive and C is positive. 
In cases where the chiral channel Ag is repulsive, or C < 0, we will explicitely set ai = 0-3 as 
the field as = (o"i — o-3)/V^ is assumed not to condense. 

In spite of the complications present in the initial formulation of this model and its 
interactions, the potential Q of Eq. (^) has a very simple structure. In particular, we 
note that it depends only on three coupling parameters. A, B, and C. The arguments in the 
logarithms of Eq. (123|) represent squares of typical excitation energies in the system. Because 
of color mixing in the diquark channel, these excitations lead in the case of two colors to a 
gap A in addition to the dynamical mass Ui. The remaining colors develop a mass 0-3, which 
is in general coupled to ui. The question of whether chiral or color symmetry breaking is 
preferred is now determined by the relative strengths of the a and A fields. This strength 
ratio depends solely on the ratio B/A or, equivalently, on the ratio of the Fierz projections 
of the original interaction onto diquark and chiral channels, B/A ~ f^/ft- The ratio C / B 
plays a secondary role, as it determines the symmetry breaking pattern only when all three 
auxiliary fields acquire non-zero mean field values. 

In the large A^ limit, the thermodynamically stable phases are determined by the saddle 
points of the potential fi, Eq. (^31) . We find that Vt has precisely one local minimum for each 
fixed value of the coupling parameters and the temperature. There are thus only second 
order transitions between single phases. The variation of B/A, C / B, and T reveals a rich 
phase structure which is illustrated in Fig. 1. Choosing the ratio B/A fixes the slope of a 
straight line in parameter space passing through the origin. Increasing the temperature T 
corresponds to moving along that line from the origin to large values of the axis coordinates. 

It is useful at this stage to repeat our objective. The only interaction that represents 
QCD is single-gluon exchange, which is a color octet, flavor diagonal, and Lorentz vector 
interaction. This interaction realizes a ratio B/A = 3/4 for Nc = 3. This correponds to 
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the dot-dashed hne in Fig. 1. All other interactions break at least one of the symmetries of 
QCD. However, a study of these deformations of QCD is useful on two respects. First, as 
discussed shortly, it reveals that the phase diagram is mostly shaped by color rather than 
chiral symmetry. Second, as discussed in Section [V^, a comparison to the phase diagram of 
microscopic models shows what aspects of the phase transition are protected by symmetry. 
We now turn to a description of the phase diagram assuming that the parameters A, B, 
and C, are chosen freely. We will later consider what ranges of values these parameters can 
actually take, given that they follow from a single primary interaction, Eq. (pl|). 

We find the usual chiral phase transition whenever B/A < Nc/2. In particular, A = and 
0"! = (73. The potential Q then has the saddle-point solutions found in the xRMM at finite 
temperature [|l^. The system breaks chiral symmetry spontaneously at low temperatures 
T < T^ with T^ = N(./{ti'^B). Below T^, we find chiral fields with a square root dependence 



in the temperature, ai = JNc/BJI — T'^/T^ |T^. Above T^, chiral symmetry is restored 
and o"! = (73 = 0. Spontaneous breaking of color symmetry occurs only for ratios B/A larger 
than Nc/2. This region of parameter space contains three phases. First, large values of 
B favor chiral symmetric phases, ai = cx-^ = 0. In this case, color is spontaneoulsy broken 
below the critical temperature Ta = 2/{n'^A'^), and is restored above. We again find a square 



root dependence for the pairing gap below Ta, A = J2/AJI — T"^ /T\. Second, keeping A 
fixed to values A < 2/(7r^T^) and decreasing B from large values, one encounters a phase 
of mixed symmetry breaking where all condensates are developed and where o"i 7^ (T3. The 
upper boundary of this phase depends on C; we have plotted three examples of this critical 
line in Fig. 1. 

Before further discussing the nature of the mixed symmetry breaking phase, we must 
distinguish cases where the chiral-Ag channel is attractive and C > from those where 
it is repulsive and C < 0. In the repulsive cases, we argued before that ai = 0-3. The 
upper boundary is then the line B = A + {Nc — 2)/7r^T^. Below that line and above the 
line B/A = Nc/2, we again find square root dependence for both chiral and diquark fields. 
Attractive cases (C > 0) require more care. Anticipating the following longer discussion 
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of those regions of parameter space which are accessible, we remark that there are strict 
hmits to the parameter ratio C / B which can be reahzed. We choose A^c = 3 for the present 
discussion. A color octet interaction gives a ratio C/ B = fy/fs = —16/3; a color singlet 
leads to C/B = f^/fi = 2/3. It is therefore clear from the definition of B and C, Eq. (pTD, 
that any combination of octet and singlet interactions which produces a positive ratio C/B 
must also satisfy C/B > 2/3. 

The pattern of symmetry breaking now depends on the interaction. Let us start by the 
case C/B = 2/3, realized by a color singlet interaction. The minimum of Q occurs for ai = 
and, since /32 = (Eq. (^l])), a^ decouples from A. The diquark field exhausts the strength 
alloted to the condensing colors, and chiral symmetry breaking acts independently on the 
third color. The upper boundary for this symmetry breaking pattern is the horizontal line 
B = 3/(7r^T^). Let us now increase C/B above 2/3. The upper boundary moves down as 
illustrated in Fig. 1 and reaches the line B = A+{Nc—2)/{tt'^T'^) in the limit C/B -^ oo. For 
any intermediate value of C/B < 2/3, the chiral fields inside the wedge of mixed symmetry 
breaking satisfy aa/ai = 1 + (3 — 2A/ B)/{C / B — 2/3). The sigma fields have again a square 
root dependence on the temperature. 

It is remarkable that the overall structure of the phase diagram depends primarily on the 
critical ratio B /A ~ Nc/2. This follows because chiral condensation involves all N^ colors, 
while quark Cooper pairs mix two colors only. 

Now comes the question of understanding whether one can really access the richness of 
the phase diagram of Fig. 1 with a random matrix model. We noted above that there are 
strict limits on the parameter ratios which can actually be realized. As we shall demonstrate 
in the following section, the upper limit on the ratio B/A is actually Nc/2. This precludes 
the exploration of phases with A 7^ 0, and thus none of the random matrix model can 
support stable diquark condensates at zero chemical. It is possible to come close to a 
stable diquark phase in the case of an axial or a tensor interaction in the color singlet 
channel. In this limiting case, the diquark condensate can always be rotated into the chiral 
condensate for colors 1 and 2, while the remaining N^ — 2 colors break chiral symmetry 
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independently. Even in this extreme case, the diquark condensate does not represent an 
independent thermodynamic phase. 

V. ALLOWED AND FORBIDDEN REGIONS OF PARAMETER SPACE 

We now wish to demonstrate that the maximum attainable ratio B /A is Nc/2. The 
quantity B/A is a weighted ratio of the Fierz coefficients of the original interaction in chiral 
and diquark channels, Eq. (l2l|) . These coefficients were obtained from projecting certain 
chiral blocks, see Eqs. ([l^) and (0). Let us ffist consider a given Lorentz channel C, 
described by ric matrices Tc {ric = 1,1,4,4,6 when C = S,P,V,A,T, respectively). Its 
Fierz coefficient in the chiral channel is /" = {nc/2)c°', where the color factor c" is c^ = 1/Nc 
for a color singlet and c*^ = 2{N^ — 1)/N^ for a color vector interaction, respectively. The 
Fierz coefficent in the diquark channel is 

ra^_. V Tr[(C7^)r^(C7^)rg] TT[{tX2)Aa (a2)Ar] .... 

where J2ii=i,nc ^^ ^ sum. over the nc members of the Lorentz channel. The ratio of traces of 
Dirac matrices can be determined from their transposition property. Noting that 

(C7')r^(C7') = -r^ when C = S,P,V, (26) 

(C7')r^(C7') = +r^ when C = AT, (27) 

we deduce that ratio to be —nc/4: when C = S,P,V and nc/4: when C = A,T. The ratio of 
traces of color matrices can be deduced from the completeness relations for the generators 
oiSU{N), 

^'-1 ^ „ / 1 \ 

H -^Q/3 Ks = 2 f 5a5 (5/37 ~ ^'^"/3 '^7'5 ) • (28) 

a=l ^ ^^c ^ 

Thus, color contributes a factor /^ ~ ~l/2 for a singlet and /^ ~ 1 + 1/Nc for a vector 
interaction, respectively. For reference, we note that single-gluon exchange produces a ratio 
of Fierz coefficients f^/fx ~ '^1^- 
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We now combine several Lorentz and color channels and try to obtain the largest ratio 
B /A. Large ratios are obviously realized by attractive quark-quark interactions. These 
channels have positive coefficients /^. They correspond to a Lorentz scalar, pseudoscalar, 
or vector with a color interaction in the A^^ — 1 channel or to a Lorentz tensor or pseudovector 
with a singlet color interaction. Their combination gives a ratio 

B ^ (Sgfo + S-f^ + 4S^f^) (iV, + I) IK + {^^-^s + ^^T\s) 1/2 

It is now clear that B /A is bounded from above as B /A < A^c/2, for Nc > 2. The upper 
bound can be reached by keeping only the last terms in numerator and denominator. This 
corresponds to either an axial or tensor interaction in a color singlet channel. 

It is instructive to consider the phase diagram of a model with the limiting value B/A = 
Nc/2. We noted earlier that the coupling between chiral fields vanishes when the interaction 
is a color singlet: now, C/B = 2/3 which implies /?2 = 0, see Eq. (J2^). The phase diagram 
follows immediately. There is only one critical temperature T^ = 2/{7i'^A) below which chiral 
and color symmetries are broken, and their fields satisfy o"| = A^ + o"^ = {2/A){l — T'^/T^). 
The symmetry breaking pattern for color is thus SU{Nc) -^ SU{2) x SU{Nc — 2). Colors 
1 and 2 exhibit a diquark condensate, while chiral symmetry operates independently on 
the third color. However, one can rotate a given color broken solution (A,(Ti) = (Ao,0) 
into (A,o"i) = (0,Ao). This rotation brings the original state to one for which ai = a^ 
and A = 0. The initial state is thus equivalent to one of broken chiral symmetry, and the 
diquark condensate does not describe a thermodynamically independent phase. We show 
in the appendix that the chiral and diquark order parameters can be described by Banks- 
Casher formulas. However, because of the rotational symmetry, all three order parameters 
are actually related to the same spectral correlator. 

We briefly return to the issue of implementing more flavor structure in the Dirac operator 
by including the flavor generators Tj. The additional flavor symmetry would contribute to 
the Fierz projections in the chiral and diquark channels. However, these additional factors 
would be those of an SU{2) symmetry, and would therefore be equal in both channels. An 
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additional flavor structure in the interaction cannot raise the ratio B/A above Nc/2. 

VI. COMPARISON WITH A MICROSCOPIC MODEL 

The phase diagram that we have obtained shares many features with those of the mi- 
croscopic model of Berges and Rajagopal 0, and other NJL models 0. We can summarize 
the three main ingredients of these models as follows. First, the four-fermion interaction 
respects the global chiral and color symmetries of QCD. Second, the coupling constants in 
the chiral and diquark channels are left as free parameters; particular choices of these param- 
eters correspond to an interaction induced either by instantons or by single-gluon exchange. 
Third, the effects of asymptotic freedom are implemented by form factors whose actual form 
is model-dependent. 

A comparison bewteen NJL-models and the present matrix models clearly shows what 
predictions on the phase diagram are protected by symmetry. To illustrate this, let us start 
from the effective potential given in Ref. [Q], and take the limit of zero chemical potential. 
With a few obvious renaming of variables, we obtain 

n = AA'' + Ba' - 2| ^|(iV, - 2){e^ + 2Thg (l + exp[-EjT])} 

+2{Ea + 2T\og (l + exp[-EA/r]) }|, (30) 

where the single particle energies are E"^ = q"^ + E^a"^, E\ = g^ + F^(cr^ + A^), and E = 
■^qcd/(9^ + ^qcd) i^ an adhoc form factor. Following Berges and Rajagopal, we do not consider 
here that the chiral condensates can split in color space, and we thus retain a single chiral 
field, 0"! = (73 = 0". Remarkably, the potential Q gives a phase diagram identical to that 
of Fig. 1. We find again a single minimum to Q for each set of parameters. This results 
in single phases separated by straight line boundaries of second order transitions. The 
general topology of the phase diagram is mostly set by color symmetry. The overall scale is 
however sensitive to form factors. Let us for instance choose the coupling constants such that 
B/A < Nc/2. For these ratios, our model has one critical line B = Nc/i'rc'^T'^) where chiral 



symmetry is restored. This line translates in the NJL model into the line B = NcJ-'{T), 
where 

^(T) = 2/^Fnanh[^] (31) 

Thus the scale l/(7r^T^) translates into T{T), which is sensitive to the particular choice of 
the form factor F . This of course implies that the actual value of the critical temperature 
is model-dependent. But the topology of the phase diagram, and in particular the line 
B = NcA/2, is protected by symmetry. 

Our approach can also be compared to the 0(N) symmetric two-dimensional model of 
Chodos et al. ||16[. This model is a generalization of the Gross-Neveu interaction to include 
pairing forces. Repeating the discussion of the previous section, one would infer a phase 
diagram similar to that of Fig. 1. Given the interaction of Ref. |T^, the critical line B = 



ANc/2 become A — B = 0. In Ref. |]T6[, A — B actually corresponds to a parameter 6 which 
determines which of the chiral {A — B > 0) ot the diquark {A — B < 0) condensates exists 
in the vacuum. The parameter 6 is further shown to be invariant under renormalization. 
This result is consistent with the fact that the topology of the phase diagram only depends 
on global symmetries. 

Although the present discussion indicates that the random matrix and NJL-type micro- 
scopic approaches share fundamental features, they differ in one essential point. NJL-studies 
generally assume from the start a potential of the form of Eq. ( ]23|) on the basis of symmetry 
considerations. The coupling contants A, B, and C, may be related to one another in specifi- 
cally motivated cases, but are in general taken as free parameters whose range of variation is 
dictated by the phenomenology to be reproduced. By contrast, random matrix models start 
out at a more microscopic level. We first selected certain types of interactions satisfying 
global symmetries, and then mimicked the dynamics by means of a random background. 
We also confined ourselves to Hermitean matrix models, in which a meaningful relation be- 
tween order parameters and spectral properties can be found. As a result of working at a 
microscopic level, the coupling parameters of the potential in Eq. (|23|) satisfy strong con- 
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straints such as B /A < Nc/2. These constraints are not deduced from phenomenology but 
are now inherited from the dynamics of the primary interactions which have been integrated 
over. Starting from a more microscopic level thus has the clear advantage of implementing 
those dynamical constraints, while still permitting the construction of the desired global 
symmetries. 

When comparing with microscopic models, it is also necessary to consider exchange ef- 
fects. Random matrix models include only direct terms. Exchange terms involve quark fields 
with unequal matrix indices k and / in Eqs. (^ and (|l^), and are therefore 1/A^ suppressed. 



To take exchange effects explicitely into account, we must modify by hand the four-fermion 
interaction produced by integration over the random background. This amounts to adding 



to the original four-fermion interaction its Fierzed transformed expression ||T5[. One way of 
proceeding is to fine tune the variances of the random background in order to reproduce 
the modified four-fermion interaction. It turns out that only scalar interactions produce ex- 
change effects that alter the chiral channel; exchange actually changes the sign of the Fierz 
coefficient in that channel. The modified interaction is repulsive for chiral condensation and 
may thus favor diquark condensation. This opens the possibility of exploring the region of 
the phase diagram above the line B /A > Nc/2. However, the matrix model that would 
mimic exchange is non-Hermitean, a case which is oustide the scope of this paper. 

A final comment concerns the similarities between single-gluon exchange and instanton- 
induced interactions in regard with the explored region of the phase diagram. To mimic 
single-gluon exchange we choose a vector interaction, which produces B/A = 3/4. To un- 
derstand how to mimic an instanton-induced interaction, we first notice that it conserves 
chirality (we assume here an interaction arising by scattering on single instantons). Thus, 
a genuine matrix model would combine scalar, pseudoscalar, and tensor interactions, and 
would then be very different in character from single-gluon exchange. It turns out that an 
instanton-induced interaction actually produces the same ratio B/A = 3/4 as single gluon 
exchange, provided exchange effects are taken into account [^ . This result is not surprising. 
The instanton-induced interaction is generated by integration over a classical gluon back- 
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ground which couples to quarks as a color vector. A consistent treatment must remember 
the vector nature of the interaction and must produce a ratio B /A = 3/4. Therefore, QCD 
and any of its consistent approximations should lie on the line B/A = 3/4. 

VII. CONCLUSIONS 

In this paper, we have introduced a random matrix model that can in principle admit 
both chiral and diquark condensation and have studied the competition of these two forms of 
order at zero chemical potential. We have considered interactions that are Lorentz invariant 
in the vacuum and gauge symmetric and have displayed the rich structure of the ensuing 
phase diagram in parameter space. The topology of the phase diagram is mostly governed 
by color symmetry. Its exploration requires the variation of only two ratios of Fierz coupling 
constants. We have further argued that there exist strong constraints on the values that these 
ratios can actually achieve. For interactions represented by Hermitean matrices, the ratio of 
the coupling constant in the chiral channel to that in the diquark channel is necessarily less 
than the critical value Nc/2 required to favor the formation of Cooper pairs. Thus, none of 
the present random matrix models can support stable diquark condensates. 

Our arguments are primarily based on symmetry considerations that also apply to QCD. 
This leads us to the conclusion that no mean-field treatment of QCD with two light flavors 
can support independent diquark condensates at zero chemical potential. This conclusion 
does of course respect the phenomenological evidence that the QCD vacuum spontaneously 
breaks chiral symmetry, but not color symmetry. The main message of our approach, how- 
ever, is to emphasize that the coupling constants of an effective potential are constrained 
by the underlying symmetries of the interaction they mean to represent. Determining these 
constraints in a given physical situation helps rule out certain forms of order on the basis of 
symmetry alone. 
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APPENDIX 



In this appendix, we explore further the case of an axial interaction with a color singlet 
structure. As we showed in Section 0, this interaction achieves a ratio B /A = 3/2. We 
found before that the gap equations yielded solutions with a finite diquark condensate which 
are degenerate with those of pure broken chiral symmetry. We now wish to derive Banks- 
Casher relations for the condensates. We show that the degeneracy of the solutions forces 
all condensates to be related to the same component of the Dirac spectrum. 

The two chiral and the diquark condensates can be obtained from the partition function 
by taking 

1 d log Z{T) 






lim lim 

mi— »0 N—>oo 

lim lim 

ms^o N^oo NNf dm-i 

1 d log Z{T) 



ml'^o N-Z^ 2NNf dnii 
1 d log Z{T) 



when a = 1, 2, 

r)=0,m3=0 

when a = 3, 

r;=0,r?ii=0 



i lim lim 



drj* 



(32) 
(33) 
(34) 



mi=0,m3=0 

Here, as usual, it is important to take the thermodynamic limit N ^ oo first to obtain 



physically meaningful quantities for the condensates |T^. The derivatives can be evaluated 
in a formal way. Integrating the partition function in Eq. (|T]) over the random background, 
an operation which we represent here by brackets, we obtain 

Z{T) = (expTrlogiP) 



where T> is the fiavor block matrix 

V = 



H + T + im rjP^ 

-r]*Pl -n^ + T -im 
Taking the derivatives in Eqs. (0), (^1), and (|3^) , we have 



V 



(35) 



(36) 



(V'LV'la) 
(^L^la) 



lim lim , , 

m^-^ON^ooNNf \(det2'D) 

,. ,. i I detzT' 
lim lim ( 

^3^0 Tv^oo A^A^^ \(deti'D) 

,. 1. -i I detzl) 
lim lim — ( ^-, ^rr ir 



*oo A^ \ (detiP) 



Tr [Pi (Pri' - ^22' 

/ ml=0,m3=0 



[a 



1,2) 



77=0,7713=0 



r?=0,mi=0 



(a = 3), 



(37) 
(38) 
(39) 
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where the color operators Pi = diag(l, 0, 0) and P3 = diag(0, 0, 1) project onto colors 1 and 
3 respectively. To evaluate the matrices "D^i, 1^22 ^ ^^^ "^uy '^^^ makes use of the fact 
that 7i is color diagonal, 7i = diag{H,H,H), and that H satisfies C'-f^ H^ C'j^ = H. A 
diagonalization of the fiavor block matrices gives then 

,,■!•,> -, , 1 / deti'D mi \ , , , , , 

{4aAa) = Im lim ^ ( j-r-^ .^^ '^ 2 ) « = 1> 2 , 40 

J- 1 / detiV m^ \ , , 

Wia^ia) = lim lim — ( , , — 2 / " = 3 , 41) 

ma^OAf^ooA' \ (detiP) [H + iy + mi,/ 

{^IP^i,,) = I hm lim ^ l^P%r , ,, J?, , „ \ , (42) 

^^^ ^^' r,^oiv^ooiV \(det2r') (i/ + r)2 + |r7|2/' ^ ' 

We thus see that the three order parameters are related to the same spectral properties of 
the Dirac operator. This result confirms the degeneracies of the solutions with finite diquark 
fields with those of broken chiral symmetry. 



ACKNOWLEDGEMENTS 

We thank G. Baym, G. Carter, K. Splittorff and J. J. M. Verbaarschot for very useful 
discussions. 



23 



REFERENCES 



[1] For recent discussions of the QCD phase diagram, see K. Rajagopal, |hep-ph/ 99083601 ; 



M. A. Halasz, A. D. Jackson, R. E. Shrock, and M. A. Stephanov, J. J. M. Verbaarschot, 
Phys. Rev. D58 (1998) 096007, and references therein. 

[2] J. Berges and K. Rajagopal, NucL Phys. B538 (1999) 215. 

[3] For recent reviews, see A. Ukawa, NucL Phys. A638 (1998) 339; E. Laermann, NucL 
Phys. B63 [Proc. Supp.] (1998) 114. 

[4] B. Barrois, NucL Phys. B219 (1977) 390. D. Baihn and A. Love, Phys. Rept. 107 (1984) 
325. 

[5] M. Alford, K. Rajagopal, and F. Wilczek, Phys. Lett. B422 (1998) 247. 

[6] T. M. Schwarz, S. P. Klevansky, and G. Papp, |nucl-th/9903048 . 



[7] G. W. Garter and D. Diakonov, Phys. Rev. D60 (1999) 016004. 

[8] R. Rapp, T. Schafer, E. V. Shuryak, and M. Velkovsky, Phys. Rev. Lett. 81 (1998) 53; 
|hep-ph/9904353 . 



[9] E. V. Shuryak and J. J. M. Verbaarschot, NucL Phys. A560 (1993) 306; J. J. M. 
Verbaarschot, Phys. Rev. Lett. 72 (1994) 2531; Phys. Lett. B329 (1994) 351; For a 
review, see J. J. M. Verbaarschot, Piep-ph/9902394 

[10] T. Banks and A. Gasher, NucL Phys. B169 (1980) 103. 

[11] E. Brezin and A. Zee, NucL Phys. B402 (1993) 613; A. D. Jackson, M. K. Sener, and 
J. J. M. Verbaarschot, NucL Phys. B479 (1996) 707; G. Akemann, P. H. Damgaard, 
U. Magnea, and S. Nishigaki, NucL Phys. B487 [FS] (1997) 721; E. Kanzieper and V. 
Freilikher, Philosophical Magazine B77 (1998) 1161; M.K. Sener, J.J.M. Verbaarschot, 
Phys. Rev. Lett. 81 (1998) 248; K. Splittorff, NucL Phys. B548 (1999) 613. 

[12] A. D. Jackson and J. J. M. Verbaarschot, Phys. Rev. D53 (1996) 7223. 

24 



[13] D. Diakonov and V. Petrov, in Quark Cluster Dynamics, Lecture Notes in Physics, 
edited by K. Goeke, P. Kroll, and H. Petry ( Springer- Verlag, Berlin, 1992) p. 288; A. 
Smilga and J. J. M. Verbaarschot, Phys. Rev. D51 (1995) 829. J. B. Kogut, M. A. 
Stephanov, and D. Toublan, |hep-pJi/9906346 . 



[14] D. Diakonov, H. Forkel, and M. Lutz, Phys. Lett. B373 (1996) 147. 
[15] See for example S.P. Klevansky, Rev. Mod. Phys. 64, 649 (1992). 



[16] A. Chodos, F. Cooper, W. Mao, H. Minakata, and A. Singh, |hep-ph/9909296 



25 



FIGURES 



Color broken, chiral 



symmetric phase 



C/B = 2/3 



C/B=l 

, Color and 
K — ---' 



PQ 



C/B-> 00 








Complete restoration of chiral 
and color symmetries 



/ 



X 



chiral brolgen 

/ 
phas6 



i^ 









i^^^-^ Chiral broken phase 







1 



2 2 

ATI T 



FIG. 1. The phase diagram in parameter space. The horizontal axis represents the coupUng 
constant in the diquark channel, the vertical axis represents that in the chiral channel. The solid 
curves are second-order phase transitions. The long-dashed lines correspond to the upper boundary 
of the phase where both color and chiral symmetries are broken. We show this boundary for various 
values of the ratio of coupling constants in the chiral and the chiral-Ag channels (see text). The 
dot-dashed line is the line of constant ratio B/A = 3/4 obtained with single gluon exchange and 
Nc = 3. This line does not cross any phase with a finite diquark condensate. 
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